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ABSTRACT—Two species of Beauveria were identified from isolates obtained from endemic 
Argania spinosa forests in Morocco, using both morphological characteristics and molecular 
data. Although isolates exhibited similar reproductive structures, ITS-rDNA_ based 
phylogenetic analysis grouped the Moroccan isolates into two clades with known sequences 
of either B. bassiana or B. pseudobassiana. Morphological comparison of the colonies 
distinguished the two groups of isolates, in full agreement with the ITS phylogenetic analysis. 
Beauveria pseudobassiana is recorded for the first time in Morocco. 


Key worps—Ascomycetes, entomopathogenic fungus, taxonomy 


Introduction 

Beauveria Vuill. (Cordycipitaceae, Hypocreales) is a genus of cosmopolitan 
entomopathogenic fungi found in various types of habitats and ecosystems 
(Meyling & Eilenberg 2007). Its species are known to infect a wide range of 
insects from many different orders in natural and agricultural environments 
(Roy et al. 2010). Preparations from B. bassiana have been developed and 
commercialized as biological control agents (Zimmermann 2007). Other 
Beauveria species hold great potential for the discovery of mycoinsecticidal 
effect (Thomas & Read 2007) and the production of different bioactive and 
pathogenetic compounds (Ames & Walsh 2010). Traditionally Beauveria 
species have been identified principally by conidial morphology. However, 
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the extensive overlap in morphological characters complicates and limits their 
usefulness for species identification. Furthermore, the discovery of cryptic 
species using molecular markers (Rehner & Buckley 2005, Ghikas et al. 2010, 
Rehner et al. 2011, Meyling et al. 2012) makes species identification a challenge 
within the genus (Fisher et al. 2011). Accordingly, accurate identification, 
including phylogenetic analyses, has been significantly improved using 
molecular methods (Rehner et al. 2011). 

Through ITS and EF1l-alpha sequence analyses, Rehner & Buckley (2005) 
demonstrated that the morphospecies B. bassiana was not monophyletic 
but comprised two morphologically indistinguishable clades that have been 
recently recognized as separate species (Rehner et al. 2011, Meyling et al. 
2012). Recently, robust molecular phylogenies based on multi-loci sequences 
have resolved 12 well-supported terminal lineages within Beauveria (Rehner 
et al. 2011). Similarly, recent sequence analyses have resulted in the recent 
description of four new species (Zhang et al. 2012, Chen et al. 2013, Agrawal et 
al. 2014, Robéne-Soustrade et al. 2015). 

In a previous study of entomopathogenic fungi in forests of endemic Argania 
spinosa (L.) Skeels (Sapotaceae) in Morocco, Beauveria and Metarhizium 
species were found but only a few isolates were identified as B. bassiana using 
molecular analysis of the internal transcribed spacer regions of ribosomal DNA 
(ITS-rDNA) sequences (Imoulan et al. 2011). The present study aims to adapt 
the phylogenetic species recognition concept (Rehner et al. 2011) by identifying 
the strains obtained from the Moroccan forests based on both morphological 
and molecular analyses. 


Materials & methods 


Fungal isolates, growth conditions, and morphological observation 

Beauveria isolates were recovered from soil samples collected from Argania spinosa 
forests in Morocco using a Galleria mellonella baiting method and maintained on Potato 
Dextrose Agar (PDA) slants at 4°C in complete darkness (Imoulan et al. 2011). The 
designation of isolates and their geographical origin are given in TABLE 1. To establish 
monosporic cultures, conidial suspensions of 1 x 10° conidia ml' were prepared from 
fungal cultures grown on PDA for 2 weeks and plated on PDA plates. The single colony 
propagated from single conidia was transferred into a new PDA dish and incubated at 
25°C. Dried agar culture vouchers were deposited in Herbarium Mycologium, Chinese 
Academy of Sciences, Beijing, China (HMAS). 

Microscopic measurements of conidia were taken from slide-cultures produced by 
inoculating a small amount of mycelium on a block of the appropriate nutrient agar 
overlaid by a cover slip. Images were acquired using an AxioCam MRc digital camera 
on a Zeiss AxioScope microscope. Microscopic measurements were performed with 
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TABLE 1. Beauveria and Cordyceps strains and ITS sequences used in phylogenetic 
analysis. GenBank numbers in bold are newly generated sequences. 


SPECIES 


B. pseudobassiana 


B. bassiana 


B. australis 


B. brongniartii 


B. asiatica 


B. sungii 


B. malawiensis 


B. vermiconia 
B. caledonica 


B. amorpha 


B. rudraprayagi 
B. varroae 


B. sinensis 

B. lii 

B. kipukae 

C. militaris 


STRAIN 


2706 
2629 

2634 

2640 

2649 

2650 

2653 

2659 

2660 

2661 

2687 

2689 

2645 
ARSEF1855 
ARSEF3405 (T) 
ARSEF6229 
2721 

2718 
ARSEF1040 
ARSEF751 
ARSEF1564 (T) 
ARSEF4622 
ARSEF4598 (T) 
ARSEF617 (T) 
ARSEF7058 
ARSEF7517 
ARSEF4850 (T) 
ARSEF4384 
ARSEF1685 (T) 
ARSEF7279 
ARSEF7760 
BCC17613 
ARSEF2922 (T) 
ARSEF2251 
AESEF2567 (T) 
ARSEF2641 (T) 
ARSEF4149 
MTCC8017 (T) 
ARSEF8257 (T) 
ARSEF8259 
RCEF3903 (T) 
RCEF5500 (T) 
ARSEF7032 (T) 
ARSEF5050 


LOCALITY 


Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Morocco 
Canada 
USA 
China 
Morocco 
Morocco 
Japan 
Vietnam 
Italy 
Australia 
Australia 
France 
USA 
Japan 
Korea 
China 
Korea 
Korea 
Malawi 
Australia 
Chile 
Brazil 
Scotland 
Brazil 
Australia 
India 
France 
France 
China 
China 
USA 
USA 


HOST/SUBSTRATE 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Soil 

Coleoptera: Scolytidae 
Lepidoptera: Tortricidae 
Coleoptera: Scolytidae 
Soil 

Soil 

Lepidoptera: Bombycidae 
Coleoptera: Chrysomelidae 
Lepidoptera: Arctiidae 
Orthoptera: Acridiidae 
Soil 

Coleoptera: Scarabaeidae 
Hymenoptera: Formicidae 
Coleoptera: Scarabaeidae 
Coleoptera: Cerambycidae 
Coleoptera: Scarabaeidae 
Coleoptera: Scarabaeidae 
Coleoptera: Scarabaeidae 
Coleoptera: Cerambycidae 
Soil 

Coleoptera 

Soil 

Hymenoptera: Formicidae 
Coleoptera: Scarabaeidae 
Silkworm 

Acari: Varroidae 

Acari: Varroidae 
Lepidoptera: Geometridae 
Coleoptera: Coccinellidae 
Homoptera: Delphacidae 
Lepidoptera 


GENBANK NO. 


KU364339 
KU364340 
KU364341 
KU364342 
KU364343 
KU364344 
KU364345 
KU364346 
KU364347 
KU364348 
KU364349 
KU364350 
KU364351 
HQ880796 
AY532022 

HQ880799 
KU364352 
KU364353 
AY531972 

AY532045 

HQ880761 
HQ880788 
HQ880789 
HQ880776 
HQ880773 
HQ880767 
AY531936 

AY532026 

AY531990 

HQ880813 
DQ376247 
HQ880824 
AY532012 

AY532003 

AY532006 

AY532008 

HQ880804 
JQ266173 

HQ880800 
HQ880801 
HQ270152 
JN689372 

HQ880803 
HQ880829 
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AxioVision Rel. 4.6 software (Zeiss, Welwyn Garden City, UK). Length : width ratios 
are given as Q, and mean values are indicated by L™ (length), W™ (width), and Q™ (Q). 


Genomic DNA extraction, PCR and sequencing 

Isolates recovered from single conidia were grown as mycelia for 2 weeks in 250-ml 
flasks containing 100ml of potato dextrose broth. Cultures were then shaken at 150 rpm 
on a rotary shaker at 25°C for 7d in darkness. The mycelial samples were pelleted from 
liquid cultures by centrifugation for 10 min at 3500 rpm (Eppendorf AG Centrifuge, 
Hamburg), subsequently washed twice with sterile double distilled water, and then 
lyophilized and stored at -80°C. 

For each isolate, 400—500mg of lyophilized mycelium was ground in a sterile mortar 
using a pestle. Total genomic DNA was extracted following the modified CTAB method 
described by Yao et al. (1999). The mixtures were incubated for 2 hours at 65°C, then 
extracted twice with 24: 1 chloroform : isoamyl alcohol, and centrifuged at 10,000 rpm 
at 4°C for 15 min. The supernatant was transferred to a clean tube and mixed with 1/10 
volume of sodium acetate (3 M, pH 5.2). Total genomic DNA was precipitated by adding 
2/3 (v/v) of cold isopropanol and chilled overnight at -20°C. DNA was recovered by 
centrifugation at 10,000 rpm, 4°C for 10 min and washed twice with 70% alcohol, dried 
at room temperature and resuspended in TE buffer (10Mm Tris-HCl, pH 8.0; 1mM 
EDTA, pH 8.0). The extracted DNA was stored at -20°C until use. 

Genomic DNA was used as template for PCR amplification of ITS region using 
universal primers ITS5/ITS4 (White et al. 1990). The PCR reactions were performed 
in a final volume of 50 ul containing 25ul 2x Taq PCR Master Mix (Tiangen Biotech 
Co., Ltd, China), 0.5 ul of each primer (10 uM), 1 ul of genomic DNA and 23 ul of 
RNase-Free water. The PCR reactions were performed in a GeneAmp PCR System 9700 
thermocycler (Applied Biosystems) as follows: 94°C for 5 min, followed by 35 cycles of 
94°C for 30 s, 53°C for 30 s, 72°C for 45 s and final extension step of 72°C for 8 min. PCR 
products were checked by electrophoresis on an agarose gel (1%) at 100 V in 0.5x TAE 
buffer (40 mM Tris-Acetic acid, pH 8.0, 1 mM EDTA), and sequenced with the same 
primer pairs. The sequencing was performed at Beijing Genomics Institute (Beijing, 
China) by using an Applied Biosystems 3730 Analyzer capillary sequencer (Foster 
City, CA, USA). All sequences generated in this study (KU364339-KU364353) were 
submitted to GenBank (TABLE 1). 


Phylogenetic and data analyses 

Only two (of 100) B. bassiana strains sequenced for ITS from Morocco were used 
as representatives, while all Moroccan strains identified as B. pseudobassiana were 
included in the analyses. Also included were 29 representative Beauveria sequences 
and one sequence of Cordyceps militaris (ARSEF 5050; as outgroup) retrieved from 
GenBank (TABLE 1). The sequences were aligned using ClustalW Multiple alignment 
tool (Thompson et al. 1997) and edited manually to avoid some obvious misalignment 
using BioEdit ver. 7.2.5 (Hall 1999). 

The sequence data set was analyzed for maximum parsimony (MP) and Bayesian 
inference (BI). MP analyses were conducted with PAUP* 4.0b10 (Swofford 2003) and 
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executed with 1000 replicates of heuristic search option of random sequence additions 
with branch swapping algorithm by tree bisection reconnection (TBR). Alignment gaps 
were treated as missing data, and all characters were unordered and equally weighted. 
Branch support was estimated by bootstrap analysis with 1,000 replicates (Felsenstein 
1985) executing the same search strategy as described above. Bayesian inference (BI) 
analyses were carried out in MrBayes 3.1.2 (Ronquist & Huelsenbeck 2003) using 
a general time-reversible model and gamma-distributed rate variation to account 
for rates of variation across sites, with a proportion of invariant sites. The remaining 
parameters were default values. One tree was saved every 1,000 generations from a total 
of 5,000,000 Markov chain Monte Carlo (MCMC) generations; the first 25% of the trees 
were discarded as burn-in. Clades with bootstrap values for MP 250% and Bayesian 
posterior probability >70% are labeled above the nodes in the phylogenetic tree. 


Results 


Morphological observation 

Microscopically Beauveria isolates from Morocco appeared typical to those 
described elsewhere (e.g., Humber 1997, Rehner et al. 2005). The colony 
characteristics were: white or pale-yellow mycelium closely (or not) appressed 
to agar surface, becoming farinaceous during sporulation and the reverse side 
uncoloured or yellowish white (Fics 1, 2). Conidiophores appeared like dense 
spherical clusters of subglobose to flask-shaped conidiogenous cells extending 
upward with an elongating sympodial denticulate rachis (zigzag appearance) 
giving rise to sessile, hyaline, holoblastic smooth conidia (Fics 3, 4). Two 
isolate groups were morphologically distinguished based on the colony colour 
as shown by aerial hyphae during conidiation: off-white and loosely floccose vs. 
yellowish-brown and densely farinaceous (Fics 1, 2). 


Phylogenetic analyses 

The amplification of ITS-rDNA region produced a uniform fragment size of 
560bp. Using Blast search, the sequences exhibited high sequence homology to 
those named B. bassiana or B. pseudobassiana in GenBank. 

The ITS dataset comprised 571 characters after the exclusion of ambiguous 
positions from both ends. There were 488 constant, 37 parsimony- 
uninformative, and 46 parsimony-informative characters. The heuristic search 
yielded six most-parsimonious trees with tree length (L) = 139, consistency 
index (CI) = 0.719, homoplasy index (HI) = 0.281, retention index (RI) = 0.807. 
MP and BI phylogenetic analyses produced trees with similar topologies 
that resolved most Beauveria lineages in separate terminal branches (Fie. 9). 
Independently of the algorithm used for phylogenetic analyses, B. brongniartii 
and B. australis were not well distinguished as unique clusters based on 
ITS-rDNA sequences. 
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Moroccan strains of Beauveria were grouped into two individual terminal 
clades corresponding to B. bassiana and B. pseudobassiana with MP bootstrap 
support >76% and BI posterior probabilities 299% (Fic. 9). These two groups 
corresponded well to their colony colours, off-white in the B. bassiana clade 
and yellowish-brown in the B. pseudobassiana clade. 


Taxonomy 
Imoulan et al. (2011) previously reported Beauveria bassiana from Morocco. 
Here we describe B. pseudobassiana as a new record for Morocco. 


Beauveria pseudobassiana S.A. Rehner & Humber, Mycologia 103: 1068 
(2011) FIGS 2-8 
Colonies growing fairly well on PDA at 25°C in the dark, reaching a diameter 
of 28-37 mm after 7 days incubation, and attaining a diameter of 36—42mm 
in 10 days; velutinous to cottony with an overgrowth of aerial hyphae up to 
the agar surface, white first and changing to yellowish-brown or pale yellow; 
becoming less powdery in older cultures and rarely with a dense farinaceous 
surface layer. Reverse side yellow with white margin. Vegetative hyphae septate, 
branched, hyaline, smooth-walled, 1-3 um wide. Conidiogenous cells solitary 
but usually consisting of dense lateral clusters, base subspherical to flask-shaped 
(sometimes ampulliform), 2.0-9.0 x 1.5-3 um, apex with an indeterminate 
denticulate rachis, produced laterally on aerial hyphae or from subtending 
cells. Conidia globose to subglobose, 2—4 x 1-3.5 um, Q = 1-2.5 (L™ = 2.6 um, 
W? = 2.1 um, Q™ = 1.4), rarely ellipsoid; hyaline, walls thin and smooth; 
produced from conidiogenous cells and occasionally directly from hyphal tips 
or laterally from hyphae. 


MATERIAL EXAMINED: MOROCCO, Ovypa; from soil under Argania spinosa tree, 23 
March 2009, coll. A. Imoulan, strain numbers 2629, 2634, 2640, 2645, 2649, 2650, 2653, 
2659, 2660, 2661 (HMAS 254118-254127, dried agar cultures); Tamanar, 28 March 
2009, coll. A. Imoulan, strain numbers 2687, 2689 (HMAS 254128-254129, dried agar 
cultures); Amskroud, 29 March 2009, coll. A. Imoulan, strain number 2706 (HMAS 
254131, dried agar culture). 


REMARKS: Beauveria pseudobassiana from Morocco differed from B. bassiana 
by its colony surface, which is yellowish-brown to pale yellow (Fic. 2) 
becoming less powdery when old, and by its occasional production of ellipsoid 
conidia (Fic. 7). The colony surface of B. bassiana is usually white to off-white 


Fics 1-8. Beauveria bassiana (strain 2718): 1. colony and aerial hyphae. Beauveria pseudobassiana 
(strain 2629): 2. colony and aerial hyphae; 3-8. conidiophores, conidiogenous cells, and conidia — 
note conidiogenous cell produced from subtending cells (6), conidia forming from a hyphal tip (7), 
and conidia produced directly from mycelium (8). Scale bars: 1, 2 = 10 mm; 3-8 = 5 um. 


Beauveria pseudobassiana new to Morocco ... 919 


920 ... Imoulan & al. 


(Fic. 1) becoming powdery when old, owing to the production of a large 
number of globose to subglobose conidia. 


Discussion 


The study of mycology in Morocco is less developed with few works being 
published on Moroccan fungi, particularly on entomopathogenic fungi 
(Imoulan et al. 2011, Imoulan & El Meziane 2013). The new species record 
for Morocco reported here demonstrates that B. pseudobassiana occurs in 
A. spinosa forests along with other entomopathogenic fungi such as B. bassiana, 
Metarhizium anisopliae, and Paecilomyces lilacinus identified previously by 
Imoulan et al. (2011). 

As shown in the most recent taxonomic revision based on the phenotypic 
criteria combined with molecular evidence (Rehner & Buckley 2005, Rehner 
et al. 2011), B. pseudobassiana and B. bassiana comprise a morphospecies 
group, morphologically indistinguishable yet phylogenetically only distantly 
related. Interestingly, our morphological observation nonetheless revealed that 
characters of both colony and conidial spores were sufficient to discriminate 
the two species. As these morphological distinctions have not been mentioned 
elsewhere, it is unclear whether our experience represents an unusual case or 
if it is possible to find additional characters to distinguish Beauveria species 
morphologically. It would be very useful if such characters could support 
DNA-based phylogenies in Beauveria. 


In the current study, ITS-rDNA-based phylogenetic analyses succeed in 
distinguishing B. bassiana and B. pseudobassiana among the species analyzed. 
Because Beauveria includes cryptic species that cannot be resolved solely 
by ITS analyses, multilocus-based phylogenies are becoming increasingly 
important for accurate assignment of strains to specific clades. Molecular 
phylogenies involving multiple genes of translation elongation factor-1a (TEF), 
RNA polymerase II largest subunit (RPB1), RNA polymerase II second largest 
subunit (RPB2), and the Bloc nuclear intergenic have become well-established 
for determining unique terminal lineages for all Beauveria species (Rehner 
et al. 2011, Zhang et al. 2012, Chen et al. 2013, Agrawal et al. 2014, Robéne- 
Soustrade et al. 2015). However, our ITS sequence analyses performed well in 
distinguishing 12 of the 15 species included, with B. australis and B. brongniartii 
grouping in a single unresolved clade, and B. rudraprayagi (MTCC8017) 
grouping with B. pseudobassiana (Fic. 9). Although the additional genes are 
required in resolving closed related species, ITS sequence is still powerful in 
identifying most of the known species of Beauveria as shown in Fie. 9. 


Beauveria pseudobassiana new to Morocco ... 921 


2716 
2593 
2721 
7659 ARSEF 1020 B. bassiana 
TNS 
ARSEF 1564 
ARSEF 751 
ARSEF 197 
B. caledonica 
100100 ARSEF 2251 


ARSEF 2922 __» B. vermiconia 


soos | ARSEF 2641 
ARSEF 6169 


B. amorpha 


13 Moroccan 
strains 


ata an B. pseudobassiana 


5573 
ARSEF 
ARSEF 3205 ; 
MTCC 3017 __, B. rudrapravagi 
RSE = 
$2.99 | ARSEF 1635 —_— 
7259 ARSEF 7279 
6599 RCEF 3903 —— B. sinensis 
ARSEF 7760 ae 
—_ ~ B. malawiensis 
BCC 17613 
ARSEF 7032 ——-» B. kipukae 
v35 ARSEF 6366 B. asata 
ARSEF 285 
zas | ARSEF S38 | B australis 
su ARSEF 2622 
ARSEF 617 
ARSEF 7053 B. brongniantii 

ARSEF 7517 

seu) ANS O57 |B varroas 

ARSEF 3259 
RCEF 5500 ———» B. lii 
ARSEFS050 Cordyceps militaris 
— image 


FIG. 9. Phylogenetic tree of Beauveria based on Maximum parsimony and Bayesian inference of the 
ITS-rDNA sequences. Bootstrap values 250% and posterior probabilities >70% are labeled above 
branches and separated by /. Terminal clades are labeled according to ARSEF accession numbers 
of individual isolates reported in Rehner et al. (2011) except RCEF5500 and RCEF3903 (Zhang et 
al. 2012, Chen et al. 2013). 
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